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The analytical expressions for both diagonal and off-diagonal ac and dc conductivities of graphene 
placed in an external magnetic field are derived. These conductivities exhibit rather unusual behav- 
ior as functions of frequency, chemical potential and applied field which is caused by the fact that 
the quasiparticle excitations in graphene are Dirac-like. One of the most striking effects observed 
in graphene is the odd integer quantum Hall effect. We argue that it is caused by the anomalous 
properties of the Dirac quasiparticles from the lowest Landau level. Other quantities such as Hall 
angle and Nernst signal also exhibit rather unusual behavior, in particular when there is an excitonic 
gap in the spectrum of the Dirac quasiparticle excitations. 

PACS numbers: 73.43.Cd,71.70.Di,81.05.Uw 



I. INTRODUCTION 

There is significant progress in fabrication of free-standing monocrystalline graphite films with thickness down to 
a single atomic layer 0, and the relatively thick (thicker than 3 monolayers) graphite films 0, 0, ! Q are now 
widely produced. The new one layer material, called graphene, possesses truly remarkable properties both from a 
technological and theoretical point of view. Graphene is a promising candidate for applications in future micro- and 
nanoelectronics due to its excellent mechanical characteristics, scalability to the nanometer sizes, and the ability to 
sustain huge (> 10 8 A/cm 2 ) electric currents. By using the electric field effect 0, 1 11 S 0, H, it is possible to 
change the carrier concentration in samples by tens times and even to change the carrier type from electron to hole 
when the sign of applied gate voltage is reversed. Another interest in graphene is related to the fact that it represents 
a building block for the other forms of carbon, viz. graphite is a stack of graphene layers, carbon nanotubes are 
wrapped graphene layers, while fullerenes can be created from graphene by introducing topological defects. 

On the theoretical side, the conduction and valence bands in graphene touch upon each other at isolated points 
in the Brillouin zone and this results in the linear, Dirac-like (up to energies of the order of 1000 K) spectrum of 
quasiparticle excitations which makes graphene a unique truly two-dimensional "relativistic" electronic system. The 
thinnest graphite films can be described by a low-energy (2+1) dimensional effective massless Dirac theory A 
recent observation 0, 13 of the unconventional integer quantum Hall effect (IQHE) 

2e 2 

^ = - — (271 + 1), n = 0,l,... (1.1) 

which is expected from the analytical study |llj | based on the fundamental properties of the 2 + 1 dimensional 
Dirac theory, can be considered as the ultimate proof of the existence of the Dirac quasiparticles in this fascinating 
material. A complementary numerical investigation of the Landau level structure for a hexagonal lattice model with 
the nearest neighbor and next-nearest neighbor hoppings also led the authors of Ref. 01 to the conclusion that the 
Hall conductivity is quantized according to the rule l|l.lf) . In contrast to this behavior expected for an ideal 2D 
graphene, thicker 2 to 10 layers thick films studied in Refs. |l|,|2j,|3j exhibit instead a conventional Hall quantization 
o~xy = ~4:(e 2 /h)n. 

The Dirac quasiparticles seem to be present not only in graphene, but also in the highly oriented pyrolitic graphite 
(HOPG), single crystalline Kish graphite and the relatively thick (thicker than 3-10 monolayers) graphite films, 
where warping introduces other types of carriers The Hall effect features in HOPG graphite were observed in 
Refs. [Til [l3 fsee also the latest Ref. 0|), but the Hall conductivity quantization in these systems remains conven- 
tional 0|yiEH E3- Nevertheless, the presence of the Dirac quasiparticles can be detected using other experimental 
techniques. For example, the differences between the Dirac and Schrodinger (massive) quasiparticles may be observed 
in thermodynamic and magnetotransport measurements 0, HI, Bl HI El E3>0] • F° r instance, the phase of de Haas van 
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Alphen and Shubnikov de Haas oscillations for Dirac quasiparticles is shifted by tt 0, H, ITM Il9l l20j compared to the 
phase of non-relativistic quasiparticles. Moreover, the Dingle and temperature factors in the amplitude of oscillations 
explicitly depend on the carrier density in the case of a Dirac- like spectrum 0] . These two characteristic fea- 
tures allow one to distinguish the Dirac quasiparticles from very light, ~ 0.01m e (m e is the electron mass) particles, 
and from other carriers which are also present in graphite 1201. The Landau levels in graphite are observed in high 
magnetic field using both scanning tunneling spectroscopy |2lj and infrared spectroscopy [22] • The latter allowed 
one to observe in HOPG the cyclotron resonance modes and to establish that some of them reveal a \ B dependence 
of the cyclotron frequency which is expected for the Dirac quasiparticles. Actually, this characteristic, yB, depen- 
dence shows up in various properties related to the Dirac quasiparticles. An interesting example is the magnetization 

M VB (cf. Eq. (7.4) of Ref. 0) at zero chemical potential /i which implies that the magnetic susceptibility 

X = dM/dB oc —B~ x l 2 diverges at zero field. Although the singularity of \{B — > 0) is smoothed [see Ref. [23j . 
where another system with Dirac quasiparticles is considered] by a coupling between layers, finite temperature and/or 
chemical potential, the presence of Dirac quasiparticles results in an anomalously strong diamagnetism of graphite 
0- 

The purpose of the present work is to extend the analysis made in the previous papers |l8l . Ilii. l2*fij (see also 
Ref. [2|| devoted to a so called (i-density wave state which is also described by the same low-energy Dirac Lagrangian) 
where thermodynamic and mostly diagonal dc magnetotransport properties of graphcnc were studied. We derive 
analytical expressions both for diagonal and off-diagonal ac conductivity which in contrast to the previous papers 
include a frequency dependent impurity scattering rate. Then we concentrate mostly on the dc Hall conductivity 
giving throughout derivations of the results presented in our short paper Ref. [TH ] and considering the limiting cases 
that were not yet considered. 

The paper is organized as follows. In Sec. |H] general features of the model for a single layer of graphene are 
described. In Sec. IIIII we present the analytical expressions for both diagonal and off-diagonal ac conductivities, 
including dc limits of these expressions (the details of calculation are given in Appendix lAl . Then in Sec. II VI we 
consider the dc Hall conductivity, the Hall angle and the Nernst signal are studied in Sec.0 In particular in Sec. lVBl 
we discuss a possibility of detecting a gap that may exist in the spectrum of the quasiparticle excitations of graphene. 
In Conclusions, Sec. lVII we give a concise summary of the obtained results. The extra technical details concerning the 
Hall conductivity in the clean limit are given in Appendix [Bj The equation for chemical potential fj, is considered in 
Appendix [C] and the solution of the Dirac equation in the symmetric gauge is presented in Appendix [DJ 



II. MODEL 

As discussed, for example, in Refs. 0,E1,II3> we start directly from the conventional QED2+1 Lagrangian density 



C=Y< *-(^ r ) [h°(m-i^)+ivF7 1 (fid x +dA^ + iv F j 2 (hd y + i^A c y ^ - A~\y a (t,r), (2.1) 
a=±i 

where f ff = ("0icr(i, r),ip2a(t, r)) is the four-component Dirac spinor combined from two spinors ipi a , ip2a [correspond- 
ing to K and K' points of the Fermi surface, respectively] that describe the Bloch states residing on the two different 
sublattices of the biparticle hexagonal lattice of the graphene sheet. In Eq. I|2.1|l 7" with v = 0,1,2 are 4x47 
matrices belonging to a reducible representation in 2 + 1, for example, 7^ = 03 <£> (0-3,10-2, —i&i), — ^i"f is the 
Dirac conjugated spinor, — e < is the electron charge, vf is the Fermi velocity, and a — ±1 is the spin variable. 
More generally the number of spin components Nf can be regarded as a flavor index and Nf — 2 corresponds to the 
physical case. 

The external magnetic field B = V x A cxt is applied perpendicular to the plane along the positive z axis and the 
corresponding vector potential is taken in the symmetric gauge A cxt = (—By/2,Bx/2). The energy scale associated 
with the magnetic field expressed in the units of temperature reads 

eBv^ ^ eBrw% 1 (Jf2) = gm x 1Q _ 8 ^ (m/s)jB(T)j (2 . 2 ) 

C C B 

where vp and B are given in m/s and Tesla, respectively. In the following we set h = ks = 1, and in some 
places e = c = 1, unless stated explicitly otherwise. There is some disagreement in the literature concerning the 
precise value of vf in graphene which is related to an uncertainty in the value of the nearest-neighbor hopping t. 
For numerical calculations we assume that t ~ 2.3eV, so that vf ~ 7.4 x 10 5 m/s which leads to the relationship 
eB -> (4.85 x 10 4 K 2 )£(T). Note that the latest experiments HQ indicate that v F « (1 -r- 1.1) x 10 6 m/s. 
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Since the Lagrangian (|2.1|l originates from a nonrelativistic many-body theory, the Zeeman interaction term has to be 
explicitly included by considering spin splitting = jjL— ag/2^BB of the chemical potential [i, where ^b = eft/ (2mc) 
is the Bohr magneton and g is the Lande factor. However, for the relativistic quasiparticle spectrum with the realistic 
values of vp oc 10 6 m/s and g ~ 2 the distance between Landau levels turns out to be very large compared to the 
Zeeman splitting |19|. so that in what follows we will not consider this term and just multiply all relevant expressions 
by the above-mentioned number of flavors Nf. We note, however, that the latest measurements in high fields (up to 
45T) |53 revealed a lifting of the spin and sublattice degeneracy, so that the half integer Hall quantization changes 
to the integer one for fields B > 20T. 

While simple tight-binding calculations (see e.g. Ref. (28[) made for hexagonal lattice of a single graphene sheet 
predict that /1 = 0, the real picture is more complicated and the actual value of ji in HOPG is nonzero due to 
inter-layer hopping, finite doping, and/or disorder. Moreover, a nonzero and even tunable value of fi (including the 
change of the character of carriers, either electron or holes) is possible in the electric-field doping experiments made on 
monocrystalline graphitic films 0, 0> El El 13 ■ I n our notations /i > corresponds to electrons and, accordingly, 
to the positive gate voltage, V g . 

The Lagrangian 1(2. ljl also includes a gap A, so that for B = it describes quasiparticles with the dispersion 
E(k) = — /1 ± y/ vpk 2 + A 2 . Again, this gap is zero when non- interacting quasiparticles on the hexagonal lattice 
with nearest neighbor hop p in g are considered. However, it could open as a result of poor screening of the Coulomb 
interaction in graphite [25l l29| and/or in the presence of an external magnetic field (the phenomenon of magnetic 
catalysis) 30] . The physical meaning of this gap (or a singlet excitonic order parameter) is directly related to the 
electron density imbalance between the A and B sublattices of the bi-particle h exag onal lattice of graphene [2|j . The 
opening of such a gap was already the subject of an experimental investigation [lfij and we hope that the predictions 
made in Refs. |25j, |29j will be tested again on the new thin samples that are closer to the ideal graphene considered 
in these theoretical papers. 

In contrast to the diagonal transport coefficients, the off-diagonal transport properties are sensitive to the sign of 
the product eB. Thus for the lucidity of the presentation, we begin with the expression for the spectral function of 
Dirac fermions and perform the calculation without assuming the positiveness of the product eB. 

The Green's function of Dirac fermions described by the Lagrangian 12.1(1 in an external field reads 

S(t - t', r; r') = exp ^rA rai (r')j S(t -t',r- r'), (2.3) 

where S(t — t', r — r') is the translation invariant part of S(t — t', r — r'). Its derivation using the Schwinger proper-time 
method and decomposition over Landau level poles has been discussed in many papers (see, e.g. Refs. [2a. l30t l31j > i . 
so here we begin with the Fourier transform of S(x — y) in the Matsubara representation 

S(tu; mj k) = e-T^T f>l)» S "^"" k ] , u m = n(2m + 1)T, (2.4) 
^ (i(j m ) 2 - Ml 

where T is the temperature, 

M n = ^A 2 + 2nv 2 F \eB\/c (2.5) 
are the energies of the relativistic Landau levels and 



S n (iu) m , k) = 2(iw m 7° + A) 



/2ck 2 \ /2ck 2 

p - L *{w\)- p + L "-iw\ 



/2ck 2 



with P± = (1 ± i7 1 7 2 sgn(ei3))/2 being projectors and L"{z) the generalized Laguerre polynomials. By definition, 
L n {z) = L° n (z) and Lt^z) = 0. 

In what follows we also need the retarded and advanced Green's functions that are obtained by analytic continuation 
from positive and negative discrete frequencies, respectively, S r (uj + i0 7 k) = S(ioj m — » ui + iO, k) and S a (uj — iO, k) = 
S(iuj m — > lu — iO, k). When the frequency dependent scattering rate T(uj) is included, they acquire the form 

g(*^(o,,k) = e-*6r E(-D" tTrr ^ ^ ' (2 ' 7) 

^ (u)±iT{uj)y ~ M 2 

The scattering rate L(w) is expressed via the retarded fermion self-energy, T(u>) = — lmI] R (uj) which in general depends 
on the energy, temperature, field and the Landau levels index n. This self-energy has to be determined self-consistently 



4 



from the Schwinger-Dyson equation. The exact form of this equation actually depends on the model assumptions 
about the impurity scattering, e.g. whether the impurity scatterers are short- or long-range and in many cases this 
equation is solved numerically. Exactly this kind of consideration was made for graphene in Ref. [^2, but in our 
paper we pursue another goal which is to obtain a simple analytical expression for the Hall conductivity. Accordingly 
here we chose a different strategy. In Sec. IHII we derive general expressions for both frequency dependent <j xx {Q) and 
a xy (Q) which include an unspecified frequency dependent scattering rate T(u>). However, we make an essential for 
the analytical work simplifying assumption that T(lu) is independent of the Landau levels index n. This assumption 
is justified when point-like impurity scattering is considered |33|, The expressions obtained for a xx (Q) and <r xy (£l) 
are suitable for investigation of microwave and optical conductivities in graphene. Then in Sec. IIVI we consider the 
case of constant width T = T(uj = 0) = -ImE fl (w = 0) = l/(2r), where r is the mean free time of quasiparticles and 
treat T as a phenomenological parameter. This approximation allows one to obtain rather simple expressions for the 
Hall conductivity in the various limits. 



III. GENERAL REPRESENTATION FOR ELECTRICAL CONDUCTIVITY 



The frequency dependent electrical conductivity tensor is calculated using the Kubo formula 

^m = — % — -, (3.i) 

where Iif-{Q) is the retarded current-current correlation function obtained by analytical continuation (H^(f2) = 
Tlij(iQ rn — > + i0)) of the imaginary time expression: 



Ikj(iS^) = y[dre ia " T {T T J i (T)J J (0)), Mr) = J ' d 2 rf{r,v), n m = 2nmT. (3.2) 
o 

Here V is the volume of the system, /3 = 1/T is the inverse temperature, and j*(r, r) is the electric current density 
operator 

f for) = = -eujpj^r.rfr^r.r). (3.3) 

a 

The brackets in Eq. (|3.2(l denote the averaging in the grand canonical ensemble. Neglecting the impurity vertex 
corrections, the calculation of the conductivity reduces to the evaluation of the bubble diagram 



^ r d 2 k 

IIy(ifi m ) = -e 2 v F T 2_, / / 27r \2 tr [l l S{iu>n^)l j S(iuj n - ffl m ,k)] , 
where tr includes also the summation over flavor index and S{iuj, k) reads 



(3.4) 



S{tu n ,k)= / ■ , (3.5) 

— oo 

with the spectral function given by the discontinuity relation 

A( w ,k) = ^-[5 A (a;,k)-5 fi (c,k)] (3.6) 

for S a - r (lo,]s.) defined in Eq. (|27jl . Note that the translat ion non-invariant phase of the fermion Green's function 
(12.31) cancels out in n. 

The sum over Matsubara frequencies in Eq. (|3.4(l is easily evaluated when the fermion Green's function is written 
using the spectral representation (|3.5|l . After this is done the analytical continuation is easily performed and we 
obtain 

oo 

n ti (O + i0) = e 2 ^ / dudu/ "*^- "* {U 1 I ^tr tfA{ U , k)^', k)] , (3.7) 
J uj — u> — il — tv J (2ir) z 



5 



where nj?(w) is the Fermi distribution function n F (ui) = l/(exp((w — n)/T) + 1). The representation 1)3.7)1 is suitable 
for studying both diagonal (see Refs. 0, H^, ) an d off-diagonal conductivities. In Appendix 1X1 we generalize 
the calculations of the previous papers and obtain both diagonal ac conductivity, 



0~xx(ty = 



e 2 N f 



duj[np(uj) — np(u>')]Re 



2B 



4tt 2 

+ [Hi(--B) + Si(+B) - Sa(-fl) - 3 2 (+B)] ^ 
and off-diagonal ac conductivity 



A 2 - (w + ;r) 2 

A 2 -(c^r) 2 
2B 



(3.8) 



w',r^ r') 



(7^(0) 



S 2 (-B)- 



2 A^ / sgn(eB) 
47r 2 fi 

2B 

i 2 - (w + ir) 2 



im 



c?w< [n^w) — rip (a/)] 



[np(w) + n F {uu')] 



Si(-B) 



(S 2 (-B)-S 2 (S))^ 

A 2 - (w + ^r) 2 



A 2 - (a; + if) 2 
2B 

(E 1 (-B)-E 1 (B))yj 



«-> w , r «-> r') 
fw + zr) 2 



2B 



(uj <-> w ',r <-> r') 



(3.9) 



Here ^ is the digamma function, and we denoted B = vp\eB\/c, uj' = uj + fi, T = T(u>),T' = T(u>') and introduced 
the following short-hand notations 



E 1 (±B) = 3 l (uj,uj',r,T',±B) = 
E 2 (±B) = E 2 (lo,u',T,T',±B) = 



(uj' + iT')(ui + iT) - A 2 



[uj - uj' + i(T - T')][uj + uj' + i(T + T')} ± 2B' 

(uj' - iT')(uj + iT) - A 2 

+ i(T + T')][uj + uj' + i(T - T')} ± 2B ' 



(3.10) 



To familiarize oneself with Eqs. 1)3.8(1 and 1)3. 9JI let us firstly consider the diagonal conductivity in two simple cases. 

Influence of the Landau quantization on the optical conductivity. We postpone a more comprehensive study of the 
influence of the gap A and the form of the function T(uj) on the ac conductivity for the future publication. Here to 
illustrate the behavior of a xx (fl) described by Eq. 1)3.8)1 in Fig.^we show the results only for A — and T(uj) — const. 
One can see that for B = there is Drude peak. However, when the magnetic field is applied the s pect ral weight is 
transferred from the Drude peak to the resonance peaks in the agreement with a recent experiment [22j. 



A. dc limits of the longitudinal and Hall conductivities 

Let us firstly consider the dc limit of o~ xx (Q), 



where 



oo 

a xx (p.B,T) = a xx (n^0) = e 2 N f J du)[—n F (u))]AL (uj, B, T, A), (3.11) 



1 



Al(u,B,T,A) - ., ., . , 
n z (vpeB/cy + (2uj1 ) z 



, (uj 2 + A 2 + T 2 )(v 2 F eB/c) 2 - 2uj 2 (uj 2 - A 2 + T 2 )v 2 F \eB\/c 



x I 2 - + " (J^A^r 2 ) 2 ^ 2 (3 - 12) 

uj(uj 2 - A 2 +T 2 ) /A 2 -(lu + iT) 2 
Imtp ' 



2v 2 F \eB\/c 

and — n' F (uj) = (1/4T) cosh _2 [(w — /j,)/2T] is the derivative of the Fermi distribution. Here the scattering rate r(o;) 
remains a frequency dependent quantity. The expression 1)3.11)1 was originally derived in Refs. pEL |26| (see also |34| 
for a related derivation of the thermal conductivity) under the assumption T(u>) — const. 
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Similarly to the dc expression l|3.11|l for a xx one can take the dc limit — > in Eq. i|3.9|) and obtain 

oo 

2u;T(B + T 2 + A 2 - tu 2 ) 



e 2 N f sgn(eB) f f , B(r 2 - A 2 + c 2 ) 

<j xv (jl,B,T) = / ■{ (-n F (w)) • 



b 2 + 4^ 2 r 2 



[r 2 + (A-^) 2 ][r 2 + (A + w ) 2 ] 



— oo 

2B 



B 



A 2 - {uj + iT) 2 



(uj + a?) (A 2 — (w + ir) 2 ) , / A 2 - + zr) 2 

V 



B 2 T \ 2B 

where now T' = dT (to) / du> . The term with 71f(w) can be integrated by parts and we finally arrive at 



a xy = e 2 N f sgn(eB) J du[-n' F (u)]A H (u,B,r,A), 



(3.13) 



(3.14) 



where 



i fs(r 2 -A 2 + ^ 2 ) 2^r(s + r 2 + A 2 



2tt 2 



- 2 ) 



B 2 + Auj 2 T 2 [r 2 + (A-^) 2 ][r 2 + (A + ^) 2 ] 



2luT A + uj A — lo 
arctan arctan 



2Imlnr 



A 2 - (u + iT f 
2B 



Im 



/ g(r 2 -A 2 + ^ 2 ) a 2 - (tj + ir) s 
^ s 2 + 4w 2 r 2 + 



B 



B 



1> 



r 

(oj + iT) 2 



2B 



(3.15) 



Recall that here B = v 2 F \eB\/c. Since we are considering non-interacting quasiparticles, the temperature and [i 
dependences of the conductivities (|3 . 1 1|> and Ij3.14|l are only contained in the derivative of the Fermi distribution. The 
spectral function l|3.15|l for the Hall conductivity turns out to be more complicated than the corresponding function 
(|3.12|) for a xx , therefore it is useful to consider simple limiting cases and to establish the correspondence between our 
answer and the results obtained by previous authors. 



IV. HALL CONDUCTIVITY 

A. Classical limit \eB\v% / c Y 2 , p? . Drude-Zener formula 

We begin our consideration with the classical limit B — ► (or more exactly \eB\vp/c <C T 2 ,fi 2 ), when Landau 
quantization is not essential. Using the asymptotic expansions 

i ^=( z -l)^- z+ b^ + ^ z -^ +o (^' (41) 

^ ) - lnZ -2^-T^ + 12cb + ^ 



we obtain that for A = and B — > 

1 BY (3T A - 8T 2 uj 2 - Slu 4 )Tlu - 3(T 2 + uj 2 ) 3 arctan^ 

Ah(w,b,t,o) - ^ l + B 2/^ u 2 T 2 ) 6r3u, 2 (r 2 + ^ 2 ) 2 ■ (42) 

Accordingly for the T = conductivity we find 

i u ^^ e 2 v 2 F N f eBn 2 sgn^ 1 

^5,0) = — ( 4ej , /c)2 + vr2 . r«| M |, (4.3) 

where we restored evp / c. The diagonal conductivity l|3.11[l in the same limit reads 



ao e 2 N f \ H i\ 

°0 = — . „ , (4.4) 



l + (w c r) 2 ' ;! 4^r ' 
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where the mean- free time of quasiparticles r = 1/2T enters instead of the transport time, r tr because we ignored the 
vertex corrections, and we introduced the cyclotron frequency, 

\eB\vl \eB\ H 
c\n\ cm c vf, 

Here in the second equality lo c is written in terms of an fictitious "relativistic" mass, m c which plays the role of the 
cyclotron mass in the Lifshits-Kosevich formula 0. The definition of m c shows that the chemical potential \p,\ in 
graphene acquires also the meaning of the cyclotron mass, so that the latter is easily tunable by the gate voltage 0. 
Accordingly the Hall conductivity (|4.3I) can also be written in the Drude-Zener form 

uj c t<t sgn(eS) sgn/i 

a ^ = r+T^rp — ' (46) 

This result agrees with Eq. (23) of Ref. when we take r tr = r = 1/2T and Nf — 1. As we will see later for Nf = 2 
all our results are twice bigger than the corresponding results of Refs. |32l IssT] . 
Finally we consider the relationship (see Appendix EJ) 

9 = ^k?/"^ (4J) 

between the chemical potential p. and carrier imbalance (density) p for the relativistic quasiparticles [H is restored]. 
This relationship is in agreement with Fig. 3 d of Ref. |7jand with Ref. |8j, viz. the cyclotron mass in graphene, 
m c = {nh 2 \p\/v F ) 1 l 2 is indeed ~ \/\p\- Since experiment 0] shows also that <r ~ |p| f° r B = 0, one may conclude 
that the carrier concentration dependence T(/x) ~ _1 for T \fx\. 
For the Hall resistivity one obtains from Eqs. (|4.4|l and l|4.6|) that 

p xy = - <*? 2 =— • (4.8) 

Thus for weak magnetic field we arrive at the standard expression for the Hall coefficient Rh — l/ecp which does not 
depend on the scattering mechanism and is used for measuring the carrier density. 

B. The limit T < \p\ < y/\eB\v%/c,T 

Another interesting and analytically treatable limit is T <C \/j,\ <C \/\eB\,Y. Since for T — > the derivative 
— n' F (u>) — > 6(u> — /it), the Hall conductivity a xy is directly expressed via Eq. H3.15JI and we need only 

A H (,,B, m . 0™ (li|±P - ^ (g)) , |,| « VB,T. (4.9) 



Substituting in Eq. (|4.9(l the asymptotic expansion of ip'(z) obtained from Eq. (|4 . 1|) in the limit y/] eB\ <cTwe have 

A H faB,r,o)~-£!£L. (4.io) 

Finally restoring the evp/c factor, we obtain 

4e 2 v 2 F N f eBfi 
° xv ~ 3^ 2 cr3 ' 

In the same limit the conductivity a xx is "universal" [2^, |3^, |36 

2 



,| « J\eB\v 2 F /c^Y. (4.11) 



e 2 N f 



(4.12) 



Moreover, theoretically the universal value (e 2 /27rft)(4/7r) (or e 2 /irh per each type of carriers) is expected even in 
arbitrary fields [2j|, because the n = Landau level is field independent. (We note that Nf = 2 when Zeeman 
splitting is neglected and Nf — 1 when it is taken into account). The experiments 0, however, show a bigger value 



8 



of the conductivity per carrier type, e 2 /h. Note also that in Ref. [37| it is argued that for the long-range impurities in 
graphene the weak-localization correction makes a positive contribution to the conductivity a xx that might explain 

the mentioned difference. 

In the opposite limit, -y/jeSf >fwe find from Eq. I|4.9|l that Ah — — 2/V(7r 2 r) and, accordingly, 



2iV/e 2 /i 



/i| <T< J\eB\v 2 /c. (4.13) 



This limit is important in the strong field (Hall) regime at small carrier densities. It allows one to extract the impurities 
scattering rate T studying the dependence of the Hall conductivity as a function of the chemical potential (or carrier 
density) in the vicinity of the point where the gate voltage changes its sign. 



C. Unusual quantization of the Hall conductivity in graphene 

In this section we discuss a full derivation of Eqs. (4) - (6) from Ref. Analyzing Eq. (6) of Ref. [ll| we 

demonstrated that the quantized Hall conductivity in graphene is equal to odd multiples of 2e 2 /h. Here we recapitulate 
the arguments of Ref. 1 1 1 1 ] discussing a few interesting moments not mentioned there. 

There are two ways to derive a xy in the clean limit. The first option is to take the limit T — > directly in Eq. i|3.7|) 
as was done in Ref. [2^| and the second option is to use a general expression (|3.14() . This option is considered in 
Appendix [Bl where we show that 



e 2 iV/sgn(eS) 

47T 



, u + A a — A 

tanh — — — h tanh ■ 



2T 



2T 



2 X! f tanh 



2T 



tanh 



[j- - M„ 
2T 



(4.14) 



Note that a xy is an antisymmetric function of \x. Rearranging terms in Eq. (|4. 14|1 and using that tanh(w — \x)j2T 
1 — 2n F (u>) one can rewrite Eq. (|4.14|) in terms of the Fermi distribution 



6 /S glHe ' V(2n + 1) [n F (M n ) + n F {-M n ) - n F (M n+1 ) - n F (-M n+1 )} . 



(4.15) 



n=0 



This representation for a xy is an equivalent of Eq. (18) from Ref. ,3§j derived for an ideal two-dimensional electron 
gas 



2 00 D / 1 \ 

— \ v /„ i i \ ( nonrclx „ / nonrcKi nonrel I „ i 1 

Vxy = -2-^{n + l)[n F {u n )-n F {uj n+1 )J, u n = ^l n+ 2/' 



(4.16) 



where m is the effective mass of the carriers with a parabolic dispersion law. The difference between the positions 
of Landau levels and their degeneracy for the Dirac quasiparticles and for nonrelativistic electron gas is encoded in 
the energies M n ~ y/n and Lu^° nTcl ~ (n + 1/2) with n = 0, 1, 2, . . ., and in the different factors 2n + 1 and n + 1 in 
Eqs. H4.15(l and l|4.16|l . respectively. 

Now we rewrite Eq. (|4. 14|1 as follows |2fij 



e 2 A^/sgn(ei3)sgn^ 
2^h 



(4.17) 



with the filling factor 



sgn/x vb = — 



, fi + A u — A 

tanh h tanh ■ 



2T 



2T 



2 ^2 ( tann 



tann 



2T 



2T 



(4.18) 



Since we are considering the quantized Hall conductivity we restore Planck constant h = 2irh in Eq. (|4.17|) and in 
what follows. Taking for definiteness /i > 0, A = and using that tanh(w/2T) = sgn(w) for T — > 0, we obtain from 
Eq. l|4~T7j> that 



>xy 



e 2 N f sgn(eB) 
2ttH 



\ + 2^9{n-M n ) 



e 2 N f sgn(eB) 



1 + 2 



9 

fj, C 



2h\eB\v% 



(4.19) 
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where [x] denotes the integer part of x. The usual semi-phenomenological argumentation (see e.g. Ref. [4C|1 for 
the occurrence of the IQHE states that in the presence of disorder when the dependence fi(p) becomes a smooth 
function, while the function a xy (fi) remains step- like. Accordingly, when in Eq. H4.19|) the spin degeneracy is counted 
by choosing Nf — 2, we arrive at the Hall quantization rule The classical (IC7|) and quantum (|4.19|) Hall 

conductivities coincide only for the odd fillings, vb = 2n + 1 as shown in Fig.|2J), where for comparison the dependence 
of o xx {vb) is also plotted. As expected the minima of o~ xx also occur at vb = 2n + 1, while the peaks of a xx coincide 
with the steps of o~ xy (yB)- 

The appearance of the odd integer number in Eq. rather than simply integer fillings is caused by the fact that 
the degeneracy of the n = Landau level is only half of the degeneracy of the levels with n > (see Appendix lD|l . 
The lowest Landau level in the irreducible representation of the 2 + 1 dimensional Dirac theory is special, because 
depending on which of two inequivalent irreducible representations is used it is occupied either by the electrons 
(fermions) or holes (antifcrmions) while at n > 1 there are solutions of the Dirac equation describing both electrons 
and holes (see Eq. (1041) ). The Dirac Lagrangian l|2.1ll which embeds a pair of independent K points on graphene's 
Fermi surface, is written using a parity preserving 4x4 reducible representation of 7 matrices that contains two 
irreducible representations with different parities. Therefore, the Landau levels associated with these K points merge 
into the full spectrum of the Dirac fermions in graphene. The superposition of the two spectra corresponding to 
the two inequivalent irreducible representations clearly explains a halved degeneracy of the lowest Landau level in 
graphene, because this level can be occupied by the holes from the K point (when /i < 0) and electrons from the K' 
point (when fj, > 0). This property of the n = level does not depend on whether the gap A has a finite value or 
A = 0. On the other hand, the higher levels may contain either electrons (when /i > 0) or holes (when /i < 0) both 
from K and K' points. 

The other way to explain the origin of "strange" odd numbers is to refer to, the mentioned above, positions of the 
minima of Shubnikov de Haas oscillations of a xx . Their unusual positions are caused by the phase shift of 7r between 
the quantum magnetic oscillations for the relativistic quasiparticles (see Refs. 0, 0, |2(J) and the corresponding 
oscillations for the nonrelativistic quasiparticles. The origin of the phase shift can be traced back to the different 
quantization of the relativistic, M n ~ y/n and nonrelativistic u^° nTcl ~ (n + 1/2) Landau levels [l8l |. 

One can gain a deeper insight into this by considering the operator § = (1/2) J d?r[w(t, r), S z ^(t, r)], where [, ] is 
the commutator and the matrix 

S z =I 2 ®(a 3 /2)= (j^ 2 ff3 ° /2 ). (4.20) 

The operator S generates the rotations of spinors in the plane by an angle <fi that are described by the operator 
U(cj)) = exp(i0§). It is natural to interpret § as the pseudospin operator, because we use the spinors that are related 
to the presence of two sublattices in graphene. Each K point of graphene's Fermi surface is characterized by a two- 
component spinor and there are two inequivalent K points. There is a temptation to make a direct analogy between 
the pseudospin operator § and the spin operator in 3 + 1 dimensional Dirac theory This, however, is misleading 
because the very notion of the spin in 2 + 1 and 3 + 1 dimensions is meaningful only for a massive particle. For a 
massless particle in 3 + 1 dimension instead of spin one introduces the helicity which characterizes the projection 
of its spin on the direction of momentum. Since in 2 + 1 dimensional case one cannot make rotations around the 
direction of the quasiparticle momentum lying in the two-dimensional plane, the helicity concept for massless particles 
is meaningless in this case. Indeed, one can check that the pseudospin operator § is identically zero for free massless 
Dirac particles. 

Nevertheless, below we argue that for the massless quasiparticles in graphene in an external magnetic field the 
pseudospin acquires a new meaning closely related to the Berry's phase discussed in the different context H,H3- 
One can obtain (see Appendix lD|l for eB > and Nf — 1 that 



:v ' ' A 



(bo m bQm 4m c 0mj + j^j- ( a nm a nm + ^ nm b nra cjj m C„ m C?Jj m d nm ) 



n— 1 m— — 00 



(4.21) 



where a nm (b nm ) are the annihilation operators of fermions with energies M n given by Eq. (|2.5(l (antifermions with 
energies —M n ) for K point and c nm {d nm ) are annihilation operators of fermions (antifermions) for K' point. The 
quantum number m < n in Eq. (|4.21|) reflects the degeneracy of each level in angular momentum. Interestingly in 
the limit A — > only the n — level contributes to S, so that the notion of the pseudospin is meaningful only for the 
states from the lowest Landau level, i.e. for the zero modes. 

Let us now consider the rotation by the angle <f> = 2ir of a quasiparticle state \n = 0) = 6j m |0) from the n = level. 
Here |0) is the vacuum state. For A = one can show that U(2ir)\n = 0) = exp(i7r)|rt = 0), i.e. after the rotation by 
2ir the quasiparticle state from the lowest Landau level changes its phase by ir. On the other hand, the states from 
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the levels with n > 1 remain invariant, because for A = the operator U(4>) does not contain the operators that 
can change these states. It occurs exactly due to this Berry's phase shift for massless quasiparticles from the n = 
level in the external field, that the carriers in graphene cannot be considered as simply very light carriers with a finite 
mass. Note that in 3 + 1 dimensional Dirac theory all fermionic states acquire the phase shift by tt after the rotation 
by <p = 2tt. 

If the n = level had the same degeneracy as the higher levels, the Hall conductivity would had been quantized in 
a more conventional manner, 

4e 2 

^semicond = n = 0, 1, . . . (4.22) 

which one might expect for a two-band [the first band would corresponds to the electrons with u n = M n — fi and 
the second band, accordingly, to the holes with w„ = —M n — fi\, two- valley [corresponding to K and K' points of 
graphen's Fermi surface] semiconductor. Again in Eq. (|4.22H we assumed that e, B, /i > 

It is appropriate here to mention that although the conventional quantized Hall conductance ct™ = — {e 2 /h)n with 
n = 0, 1, ... is often derived by solving the nonrelativistic Schrodingcr equation, it was shown [4ll |42| that there are 
no relativistic corrections to this expression and the same result remains valid for a relativistic electron gas confined 
in the plane described by the (3 + 1) dimensional Dirac equation. The case of graphene is different, because the Dirac 
theory is used to describe an effective theory of non-relativistic quasiparticles with a linear dispersion. Although 
one usually associates the Dirac- like description of graphene with a linear dispersion of quasiparticle excitations, 
E(k) = — /z ± «F|k| for B = 0, the quantization i|l.l|) survives even when there is a nonzero gap A. 

Another important feature of graphene is that its Dirac-like description is based on the 4x4 reducible, parity 
preserving representation of 7 matrices which allows one to include two inequivalent K points of graphene's Fermi 
surface. From a theoretical point of view one may also choose a separate K point and consider the role of the parity 
breaking terms (see e.g. Refs. 0,0])- This approach is closely related to early unsuccessful 0] attempts to explain 
the IQHE using the chiral anomaly. However, in the case of the parity anomaly the Hall effect occurs even in zero 
magnetic field and in the absence of Landau levels [43. On the contrary, the Dirac-like description of graphene 
preserves parity, so that the Hall conductivity is always absent in zero magnetic field. 



D. Illustrations of analytical results for the Hall conductivity 

Fig. [3 is plotted to illustrate the odd integer Hall quantization ijl.lfl and it is computed on the base of Eqs. (|3.14() 
and I3.15f) . Since the representation (|3.14|) . (|3.15l) is derived under the approximations discussed in Sees. ITU and 
IIIII this consideration ignores the presence of localized states. Thus for finite scattering rate T this representation 
for a xy (B, r, A) provides only an approximate description of the Hall quantization, because T remains nonzero even 
between Landau levels. Nevertheless, considering the relative simplicity and analytical character of the two expressions 
(|3.11|l for a xx and l|3.14|l for cr xy , overall they give an amazingly good description of the conductivities. In particular, 
one can see o~ xx is indeed very small in the plateaux regions of a xy . 

The IQHE can be obtained by varying either carrier concentration as done in Fig. [5] or the value of the applied 
magnetic field. The latter possibility is shown in Fig. |3 where both conductivities (Fig. a) and resistivities (Fig. [31 
b) are plotted at T = 3K. To observe the quantum Hall effect in conventional semiconductors one should go down to 
the temperatures lower than 10K, while in graphene it can be observed up to 100K. This wider range of temperatures 
is related to the fact that for the Dirac quasiparticles the distance between Landau levels is much larger than for 
the nonrelativistic quasiparticles in the same applied field, so that some signatures of Shubnikov de Haas oscillations 
remain notable even at room temperature ■ 

Although the definition of the filling factor [3^| does not depend on whether relativistic or nonrelativistic quasipar- 
ticles are considered, the relationship (|4.7|l between the chemical potential and carrier imbalance shows that the small 
filling factors become accessible in relatively small, compared to conventional semiconductors, fields. These features 
make the IQHE in graphene very promising for fundamental research and possible applications. 

Finally one can observe that a xy = for B — > 0. This illustrates the point mentioned in Sec. II V CI that the IQHE 
in graphene is conventional in the sense that its explanation does not rely on any kind of parity breaking anomaly 



E. Magnetic catalysis and its observation in the Hall conductivity 



The flat, field independent behavior of a xx (B) and a xy (B) for B > 1.1T seen in Fig. corresponds to the regime 
where only the lowest Landau level is filled, because it always stays below the chemical potential (except for /1 = 0). 



11 



As we already mentioned in Sec. [HJ it was predicted in Refs. |25t |30( that for Dirac fermions in 2 + 1 there is a 
phenomenon called magnetic catalysis. It is expected that above a critical field B c which is a function of fi and T, 
a gap A should open in the spectrum of the Dirac fermions. Since the conditions for the magnetic catalysis are 
the most favorable at low carrier concentrations (or /x 0), in Fig. 0] we present the behavior of a xy (B, A(B)) and 
&xx{B, A(-B)) for \/j,\ smaller than was used to plot Fig. EI For comparison we plot these conductivities both for A = 
and for the phcnomcnological gap dependence 



A{B) = Cy /B - B C 9{B - B c 



(4.23) 



where c is some constant. As stated above, B c — B c ([jl, T), but for illustrative purposes we simply choose an arbitrary 
value B c = 0.3T. As one can clearly see, the opening of the gap A causes the decrease of \a xy (B)\ from the last 
plateau value 2e 2 /h. This tendency can also be understood from Eq. (|4.14|) . In the strong field limit the sum over 
M n with n > does not contribute and we obtain (for eB > and [i > 0) 



e 2 N 



f 



>xy 



4ttH 



tanh — + tanh r 



2T 



2T 



e 2 N 



2ttH 



t-9 (n — A(B)), 



T — > 0. 



(4.24) 



Eq. H4.24f) shows that when A(B) > fj, the magnetic catalysis leads to a formation of a new insulating phase. On 
the other hand, we observe in Fig. 01(a) that cr xx (B) increases as the gap opens, while the increase of the diagonal 
resistivity (Fig. 01(b)) also indicates the system goes towards an insulating phase. 

In our consideration we assumed that the opening of the gap A does not affect the chemical potential /x. We will 
come back to this important issue in Sec. IV Bl 



V. HALL ANGLE AND NERNST COEFFICIENT 

The approach presented allows one to calculate other transport coefficients such as thermal conductivity Kij (B, T, A) 
[see Refs. [l9l l2o. l34| . where the diagonal thermal conductivity is studied] and Peltier (thermoelectric) conductivity 
f3ij(B, r, A) tensors. The calculation of the off-diagonal coefficients involves some subtleties, because the conventional 
Kubo expressions have to be altered [4ll l47l lisl ] to reflect the role of magnetization on the electronic thermal transport 
in applied field. However in the low-temperature limit one may rely on the Sommerfeld expansion and express the 
thermoelectric tensor through the conductivity tensor 



7T 2 Tda l3 
l l3 ~ 3 e cV ■ 

The Nernst signal measured in the absence of electric current is expressed in terms of and /3y- as 



e y (T,B) 



Ey 
V X T 



&x x fix y y Px 

a 2 + a 2 

^ xx 1 ^ xy 



(5.1) 



(5.2) 



Since in the low-temperature limit Eq. (|5.1|l is valid, the Nernst signal (|5.2J) can be found differentiating the Hall 
angle 



Oh = arctan 



via the relation 



e y (T,B) 



7T 2 T dO H 
3 e d/j, 



(5.3) 



(5.4) 



Here based on the results derived in the previous sections, we will study the behavior of the Hall angle and Nernst 
coefficient. 



A. Drude-Zener \eB\vp/c « T 2 « /i 2 and T < < ^\eB\v 2 F /c < T limits 
In the classical limit (see Sec. II V A} a xy and a xx are given by Eqs. (|4.6() and (|4.4() . respectively. Hence, 



Qh — tan Oh = —u c t sgn(ei?) sgnzx 



v 2 F eB 
2cTfi ' 



(5.5) 
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The Hall angle (|4.6(l diverges at n — > 0, but one should consider only finite values of /i, because it is derived under the 
assumption \eB\v 2 F /c <C T 2 <C /.t 2 . Accordingly when /x is small, Eq. (|4.(il) is valid only for very low fields B. When 
the character of the carriers changes from hole-like (fi < 0) to electron-like (/i > 0), the Hall angle also changes from 
positive to negative. The behavior of the Nernst signal 



7T 2 k%Tv 2 F B 

3 My" 



(5.6) 



mirrors the divergence of the Hall angle at fi — > 0, but e y < irrespective of the sign o f a , [Boltzmann constant ks 
is restored in Eq. I|5.6[l .] Eqs. (|5.5|l and l|5.6[l agree with the results obtained in Ref. [43 using Boltzmann theory. 
This is not surprising, because this is exactly the limit described by Drude-Zener theory. We note that although in 
Ref. [4?]j | the so called d- density-wave state is considered, a direct comparison with our case is possible, because the 
effective low-energy Dirac theory turns out to be the same |26j . 

Since in the regime T <C \fj,\ <C \J\eB\vpj c <C T the conductivities a xx and a xy are given by Eqs. (|4.12|) and l|4.11|l . 
respectively, for the Hall angle one obtains 



_ 4v%eB» 

° H --^cT^- (5J) 



Hence the Nernst signal is positive and given by 



4^ 2 k%Tv F B 

-9 J3— • V*) 



The Nernst signal Ij5.8|l can be very large in clean system because e y ~ r 3 This regime is accessible due to 

the fact that fj, « 0, i.e. one may consider the large and positive Nernst signal as another fingerprint of the Dirac 
quasiparticles. We mention that since the Dirac quasiparticles emerge in the scenarios with unconventional charge 
density waves (UCDW) (2^, H(j , in Ref. 50] the large and positive e y is regarded as a hallmark of UCDW. 



B. Illustrations of analytical results and detection of the gap A from the Hall angle measurements 

In Fig. [S] and El we present the dependence of the Hall angle (|5.3f) on fj, for two different values of the field B. The 
case of small B shown in Fig. agrees with the analytical expressions discussed above. Indeed when decreases, 
there is an increase of ®h [cf. Ea. l|5.5|) ] followed by the regime where Qh crosses zero [cf. Eq. I|5.7ll ]. 

In Fig. [7| and[S] we show the behavior of the Nernst signal. When is large and B is small, e y < and rather 
small. However, when decreases, e y becomes positive and large (~ lOOpV/K) in accord with Eq. (|5.8|l . When the 
field B increases, the value e y (fi — 0) becomes even larger and for finite [i there are oscillations of e y (p). 

Analyzing Figs. one may discover another interesting property, viz. in the presence of a nonzero gap A the 

dependence &h(^, A) in Figs. El near \i 1=3 is not so steep as compared to the A = case. Accordingly this 
feature is reflected in the Nernst signal, so that a finite A also shows up in the dependence e !/ (/i, A) as a dip. Thus we 
suggest that careful study of the Hall angle Oh(h, A) and Nernst signal e y (/z, A) may help to establish the presence 
of a nonzero gap A in graphene. 

It is important, however, to stress that the gap is detectable only if its opening does not change the chemical 
potential \x. The situation is exactly the same as in Ref. [26^. where we considered the possibility of detecting the gap 
using precise measurements of the period of the quantum magnetic oscillations (de Haas van Alphen or Shubnikov de 
Haas). In the clean system with a fixed carrier density p, the chemical potential fi is given by the number equation 
(IC6I) . The opening of the gap A results in the adjustment of fi, so that the gap cannot be detected neither from the 
period of the oscillations |26j | nor from the Hall angle. Nevertheless, the observation of the quantum Hall effect 0, 
shows that there is localization and the chemical potential remains fixed making the gap detection possible. 

Finally we stress that here we neglected the dependence of T(/i), while a simple argument given at the end of 
Sec. IIV A1 shows that this carrier concentration dependence is rather important. 



VI. CONCLUSIONS 



In this paper we have studied the dc Hall conductivity, Hall angle and Nernst signal in a planar system with 
relativistic Dirac-like spectrum of quasiparticle excitations. We also presented the results for the diagonal optical 
conductivity in the external magnetic field perpendicular to the plane. 
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Our main results can be summarized as follows. 

(1) We have obtained analytical expressions for the diagonal optical conductivity a xx (£l, /it, B, T, A) [Eq. i|3.8p ] and 
the off-diagonal optical conductivity cr xy (Q, /i, B, T, A) [Eq. I|3.9[l]. 

(2) We have derived the analytical expression (|3.14|) with the kernel (|3.15|) for the dc Hall conductivity 
a xy (^i, B, T, A) which includes an arbitrary impurity scattering rate T(lo) that was assumed to be independent of 
the Landau level index n. 

(3) We have shown that in the classical limit our expression for dc Hall conductivity Eq. (|3.14|) reduces to the 
conventional Drude-Zener formula (|4.(il) . 

(4) The direct comparison of the expression l|4.16|l derived in Ref. [38| for Hall conductivity in a two-dimensional 
electron gas with the corresponding representation 14.15fl for graphene allows one to understand the origin of the odd 
integer Hall quantization Eq. (|l.lf> in graphene in terms of the difference between the energies and degeneracies of 
the Landau levels in these systems. 

(5) In Sec. IIV O we presented the arguments (using the second quantization formalism) that the nontrivial Berry's 
phase in graphene is associated with the anomalous properties of the zero modes or the quasiparticles from the lowest 
Landau level. 

(6) We have investigated the behavior of the Hall angle and the Nernst signal showing that for // « there is an 
interesting regime [see Eq. <|5.8|l ] where the Nernst signal is strong and positive. 

(7) On the basis of the results obtained, we have discussed in Sees. HVEl and lVBl the possibility of detecting a gap 
A that may open in the spectrum of the Dirac-like quasiparticle excitations of graphene due to a nontrivial interaction 
between them. 

All our results are derived for noninteracting quasiparticles treating the impurity scattering rate L(w) as a phe- 
nomenological parameter and without considering the interaction with impurities that would demand solving an 
equation for L(w, n) (see Ref. |32^). Accordingly we did not consider the problems related to localization and a full ex- 
planation of the IQHE (see Refs. ]40H5l|). These problems by themselves acquire a new depth and deserve a separate 
study. For exam ple, it is pointed out in Ref. Q that localization effects in graphene are suppressed. Interestingly, it 
is shown in Ref. |37| that for long-range impurities the weak-localization correction makes a positive contribution to 
the conductivity <j xx . This anti-localization property is related in Ref. |37j | to the Berry's phase of the Dirac fermions. 
Definitely, such effects were not considered in the present work. Nevertheless we hope that the approach presented 
here allows one to explain in the most transparent way the difference between the Dirac quasiparticles in graphene 
and nonrelativistic quasiparticles in conventional semiconductors when these systems are placed in a magnetic field. 
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The most efficient way of calculating <Jij(Q) is to use again the spectral representation Eq. (|3.5|) (for /i = 0). This 
allows one to eliminate one of the integrations over frequency in Eq. (|3.7() , so that for the real part of the conductivity 
(|3.1|) we obtain 
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APPENDIX A: CALCULATION OF ^-(fi) 



oo 




tr { [np(uj) — uf(oj + n)]Ys R {uj + n, k)^ j s A (uj, k) 



(Al) 



— OO 



n F (cu)-/ l S R (u + Q,k)^S R (uj,k)+n F (Lu + n) 1 l S A (cu + Q, k)^S A (u, k)} . 
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Using the expressions Eq. (|2.7|l for the advanced and retarded Green's functions we obtain 

oo , 

d 2 k _. 



(27T) 



^ (-l) n+m tr<^ [n F (w)-n F (w + n)] 



7 ^(c + ft) 7 ^M 



[(w + n + iT{uj + ft)) 2 - Af 2 ][( W - - 
7 ^(cJ + ft) 7 ^H 



— n F (uj) 



-fs*{u> + tty?s*{u>) 



[(uj + ft + iT(w + ft)) 2 - Af 2 ][(w + zr(w)) 2 - Ml 



+ n F (uj + Q) 



[(w + ft - iT(uj + ft)) 2 - M 2 ] [(w - t r(w)) 2 - Ml 



(A2) 

where S^' A ^ are the numerators of Eq. I|2.7(l obtained from Eq. (|2.6() for S n (iui n ,li) via the rule UJ, K = 

S n (iu> n — > uj ± ir(w)). This allows to include the frequency dependent impurity scattering rate, T(u>). The traces in 
Eq. I|A2J) are easily evaluated 



tr 



YS^iu'WS^iu)] = -8JV)[(w' ± tT')(w ± ir) - A 2 ] [% (L n (x)L m _i(a;) + L n . x (x)L m {x)) 



+ie tJ sgn(eB) (Z^aOZm-ifa) - L n -i(x)L m (x))] ~ 64N f (2k t k 3 - S lJ k 2 )L 1 n _ 1 (x)L 1 m ^ 1 (x), 



(A3) 



where ey is antisymmetric tensor (ei2 = 1) and the argument of the Laguerre polynomials is x — 2ck 2 /\eB\. Inte- 
grating over momenta we obtain 

e 2 v 2 F \eB\N f 



(Tij (ft) 



•>-Vft Re ^ (~ 1 )" +m+ M*y(^n,m-i +^n-i,m) + ie y sgn(eS)((5„ !m _i - (5„_i^„)] 
r[n F (w) -n F (w')][(w' + ir')(^-ir) - A 2 ] n F (w)[(w' + ir')(w + ?;r) - A 2 ] 

w I [(^ + zr') 2 -Af 2 ][(c-zr) 2 -M^; " — 

ti f (uj')[(uj' — iT')(uj — iT) — A 2 ] 



[(c' + i r') 2 -M 2 ][(c + i r) 2 -Af 2 l ; 



(A4) 



[(c'-iT'^-M^^-ir^-M^] 



The Kronnecker's delta symbols appeared in Eq. (|A4I) are due to the the orthogonality relation for the Laguerre 
polynomials, 



(A5) 



dxe x L (x)L (x) — : 5 mn , 

n! 

show that only transitions between neighboring Landau levels contribute in the conductivity. The last term of Eq. (|A3 
vanished after the angular integration. The real part of Eq. 1A4JI reads 



(Tij (ft) 



e 2 v 2 F \eB\N f 

27T 2 cft 



n,m— 

eij sga(eB)(6 n>m -i 



OO 

CO „ 

E (-l)" +m+1 / dto {<%(<Wi + S n - ltm )[n F (u) - n F (w')]Re[ILl m (uj,uj') - LT 2 ^, w')] 



where we introduced 



n n,m(^^') = TTT 



*n-l,m)NH ~ nf(w')]Imn,\ m (w,w') - [n F (w) + n F (c/)]ImII 2 m (w,t/)]} , 
(u/ + zL')(lj — ir) — A 2 2 (o/ + iT')(u) + iT) — A 2 



(A6) 
(A7) 



[( W '+ff') 2 -M 2 ][(a;-ir) 2 -M 2 ] ! ' ' [(w' + iT') 2 - Af 2 ] [(uj + iT) 2 — M^ 

When we derived Eq. I|A6(I we used the fact that the real part of II^' 2 ^ does not alter when the simultaneous re- 
placements iT — » — iT and iT' — ► — iT' are made, while its imaginary part reverses sign. These features of cry (ft) are 
also used below when Eqs. (|3.8I) and (|3.9() are written in the symmetric form. The sum over m in Eq. i|A6(l is easily 
calculated using Kronnecker delta's, 



(ft) 



e 2 v 2 F \eB\N f 

27T 2 cft 



OO 

oo „ 



n 2 „+i,Jw,w')] - ey sgn(aB) [[ti f (uj) - 7i F (u/)]Im[n^ n+1 (u;,u/) - TT) l+l n (u),u)')\ 
[n F (uj) +n F ( W ')]Im[n^ +1 ( W , W ') - n 2 +M >, a/)]] } . 



(A8) 
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Further summation over n can be performed expanding II ' in terms of the partial fractions, 



n n,m( w V) 



(lu' + iT')(u) - iT) - A 2 



[w - w' - i(T + T')} [lu + to' - i(T - T')} + 2B(n - m) 

(m' + iT')(u) + iT) - A 2 

[lu-lu' + i(T- T')} [u + lu' + i(T + T')} + 2B(n - m) 



1 



1 



(lu' + iT') 2 - Ml [u> - iT) 2 - Ml 
1 1 



(lu' + iT') 2 - M 2 (w + iT) 2 - Ml 



(A9) 



where for brevity we introduced the notation B = v 2 F \eB\jc. The resulting sums are expressed via the digamma 
function by means of the formula 



E 

71=0 



.4 



B 



C 



D 



n+a n+b n+c n+d 



\Ail>[a) + B^(b) + C^j(c) + D^P(d% 



(A10) 



where for convergence A + B + C + D = 0, so that we arrive at the final expressions for the conductivities l|3.8() 
and l|3.9l) . In the limit of vanishing magnetic field the Hall conductivity becomes zero, while for the longitudinal 
conductivity we obtain 



e 2 N f f 

— oo 

(u + iT)(uj' - iT') - A 2 



(lu + iT)(lo' + iT') - A 2 



[lu-lu' + i{T + T')} [lu + lu' + i(T - V 
Calculating the real part we can represent the last expression in the form 



lu-lu' + i{T- T')} [lu + lu' + i(T + T')} 
ln[A 2 - (lu + iT) 2 } + (lu <-> lu' , T <-> T') 



&xx(ty 



e 2 N f 
2tt 2 



du) 



Hf(ijj) — uf(lu + f2) 



(All) 



(A12) 



where 



1 



D(lu,lu') 



In 



+6 arctan 



(A 2 + T 2 - ^ /2 ) 2 + 4m' 2 T" 
(A 2 +T 2 -lu 2 ) 2 +4lu 2 T 2 
2luT 



A 2 +T 2 



c arctan 



2lu'T' 



LU* 



A 2 + T' 2 



and 



D(lu, lu') = [(T - T') 2 + (lu- lu') 2 } [(T + T') 2 + (lu - lu') 2 } [(T - T') 2 + (lu + lu') 2 } [(T + T') 2 + (lu + lu') 2 } , 
a =TT' [(lu' 2 -lu 2 + T' 2 - T 2 )[(lu' 2 + lu 2 + T' 2 + T 2 ) 2 + A(lu' 2 lu 2 - T' 2 T 2 )} - 8A 2 lulu'(lu' 2 - lu 2 - T' 2 + T 2 )] 
b =2 [luT'(lu 12 +lu 2 + T' 2 + T 2 )[(lu' 2 - lu 2 ) 2 + (T' 2 - T 2 ) 2 + 2(lu' 2 + lu 2 )(T' 2 + T 2 ) - 8lu' 2 T 2 } 

-2A 2 cu'T'[(cu' 2 - lu 2 ) 2 + (T' 2 - T 2 ) 2 + 2(cu' 2 + lu 2 )(T' 2 + T 2 ) - 8cu 2 T 2 }] , 
c =2 [lu'T(lu' 2 +lu 2 + T' 2 + T 2 )[(lu' 2 - lu 2 ) 2 + (T' 2 - T 2 ) 2 + 2(tu' 2 + lu 2 )(T' 2 + T 2 ) - 8lu 2 T' 2 } 

-2A 2 cuT[(lu' 2 - lu 2 ) 2 + (T' 2 - T 2 ) 2 + 2(lu' 2 + cu 2 )(T' 2 + T 2 ) - 8lu' 2 T' 2 }] . 



For A = these expressions are simplified: 

a = rrv 2 - u 2 + r' 2 - r 2 )[(^' 2 + lu 2 + r' 2 + r 2 ) 2 + a(lu' 2 lu 2 - r' 2 r 2 )], 

b = 2 [luT'(lu' 2 +lu 2 + T' 2 + T 2 )[(lu' 2 - lu 2 ) 2 + (T' 2 - T 2 ) 2 + 2(lu' 2 + lu 2 )(T' 2 + T 2 ) - 8lu' 2 T 2 }] , 
c = 2 [lu'T(lu' 2 +lu 2 + T' 2 + T 2 )[(cu' 2 - lu 2 ) 2 + (T' 2 - T 2 ) 2 + 2(lu' 2 + lu 2 )(T' 2 + T 2 ) - 8cu 2 T' 2 }] . 



(A13) 



(A14) 



(A15) 



The behavior of a xx (Q) in this limit was studied in Ref. |52| for graphene and in Ref. |53| for a d-wave superconductor. 
For B = the latter is rather similar to graphene at n = 0. 
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APPENDIX B: DERIVATION OF THE HALL COND UCTIVITY IN THE CLEAN LIMIT FROM 

Eq. I|3.14fl 

Eq. (|4.14|) can be obtained directly from Eq. (|3.14(l . Indeed taking the limit r — ► 0, we find from Eq. (pH5|l that 



2tt 2 



arctan ■ 



2ujT 



- 2Imlnr 



A 2 + T 2 - u? - HujY 
B 



r -► o. 



(Bl) 



Recall that here B is a short-hand notation for v F \eB\/c. We begin with the expression 



ImlnT 
= sgn wlm In T 
Now we use the relationship 



A 2 +T 2 - w 2 - 2iLoT\ 
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2B 
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(B2) 



T(z)T(-z) = 



z sin(— irz) 



(B3) 



to rewrite the last expression in the form 
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(B4) 



Hence 



Ah(u, B, T, A) = | „ 4 arctan ■ 
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Now using the formula 
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O ^ 

Imlnsin(a + ib) = — - a - V" - sin(2ka)e~ 2kb , b > 0, 

fe=i 



(B6) 



we arrive at the following representation 



sgncj 

Ah(v, B, T, A) = < arctan - 
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(B7) 



Finally by means of the identity (4.6) of Ref. we obtain in the limit r — > that 



A H (w, B, T, A) = -— sgn wQ(u? - A 2 ) 



2tt 



fc=l 



sin 2itk- 
ttA; V 2B 



-— sgnw 

Z7T 



6»(w 2 - A 2 ) + 2 J] 6»(w 2 - A 2 - 2Bn) 



(B8) 



Inserting the last expression in Eq. I|3.14|l and integrating over u) we finally arrive at Eq. I|4.14[) . 
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APPENDIX C: THE EQUATION FOR CHEMICAL POTENTIAL 



To derive the relationship for the carrier imbalance p (p = n e — rih, where n e and rih are the densities of "electrons" 
and "holes" , respectively) and the chemical potential p, we begin with the well-known expression 



/7 = tr[ 7 5(r,0)], r^O, 



(CI) 



where S(t,t) is the translation invariant part of the Green's function l|2.3(l . Taking its Fourier transform and using 
the spectral representation (|3.5(l . we arrive at 



71 — - — OO 



d 2 k f°° tr[ 7 A(q;,k)] 
(2tt) 2 J_ 00 iuj n + fi-uj' 



After evaluating the sum over Matsubara frequencies we obtain 

1 f d 2 k 



P : 



tanh - ' ^ trfy A(oj , k)]. 



(C2) 



(C3) 



Now taking into account that tr[7°A(w, k)] is an even function of u>, we may write 



1 f d 2 k 



P= 2 



(27T)" .,„ 
d 2 k 



du> 



u> + a to — (jl 

tanh tanh 



trfr^k)] 



(2^) 2 Jo 



2T 2T 
du> [iif(li — ji) — np{u) + pi)] tr[7 A(u>, k)]. 



(C4) 



Substituting the spectral function .(j|> in the first line of the last equation and integrating over momenta, we arrive 
at 



N f \eB\ ^ 



Att 2 c 



a,, 



n=0 



diu tanh 



w + p 
2T \(Lu-M n ) 2 +T 2 



+ (w — > -w) , ao = 1, a„ = 2, n > 1. 



(C5) 



The ratio a„/ao = 2 for n > 1 is related to the above-mentioned smaller degeneracy of the n — Landau level. It is 
easy to see that Eq. (|C5ll in the limit r = reduces to 



Nf\eB\ 
Aitc 



, u + A u, 
tanh — — — 1- tanh ■ 



2T 



2T 



2 E( tanhJ 27 



tanh 



p - M n 
2T 



(C6) 



while in the limit B = r = T = A = 0it reduces to Eq. (|4T7|l (see, for example, Eq. (77) in Ref. giving the 
relationship between the chemical potential and the free carrier imbalance. 

Comparing Eq. i|C6(l with Eq. I|4.14|l we finally obtain that in the ideally clean system (see also Ref. [25j) 



ecp 

"IT' 



(C7) 



The last expression seems to be paradoxical at first glance, because it corresponds to the classical expression (|4.8|l far 
beyond the validity of the classical limit (see Sec. II V All . Nevertheless this result is absolutely correct and it shows the 
consistency of our calculation. As explained in Ref. 153. Ed. (|C7I) is expected for an ideal conductor. This similarity 
between Eq. 1|C6() and Eq. 1)4. 14|) was exploited in Ref. 44], where instead of calculating the electrical conductivity a xy , 
the density l|C6(l was obtained. However, to consider the IQHE one must take into account the presence of impurities 
p0ll5ll l53|. It is believed that they lead to the localization of most of the bulk states, except in a region around the 
center of the Landau band, and act as a reservoir which almost fixes the chemical potential [4(i l5l| . It turns out that 
in this more physical situation it still makes sense to rely on Eq. 14.1411 for a xy , while Eq. (|C6J) cannot be used in the 
IQHE regime. One particular model of the equation for (i that includes the reservoir could be a generalization of the 
one discussed in Ref. |55| for nonrelativistic quasiparticles. 

The main implication of this model is that the density of the delocalized carriers in the Hall bar may oscillate as 
the field B varies. It seems the oscillations of this kind were indeed observed in the IQHE system [5(|. However, there 
is no consensus on a microscopic picture of the localization in the quantum Hall effect (see e.g. Refs. [53, lla l59j 1 
even in 2D electron gas. On the other hand, the localization of the Dirac quasiparticles appears to be quite different 
from the localization in 2D electron gas with the parabolic dispersion • This indicates that further studies of the 
localization and oscillations of the density of the delocalized carriers in graphene may be very useful. 
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APPENDIX D: SOLUTION OF THE DIRAC EQUATION IN THE SYMMETRIC GAUGE 



The Dirac equation in the problem of a relativistic fermion in a constant magnetic field B takes the following form 
in 2 + 1 dimensions: 



i^ hd t + ivp^f 1 (hd x 



+ i 



■ WfJ 2 



(tidy 



1&(t,r) =0, 



(Dl) 



where the vector potential A ext = (— By/2, Bx/2), so that the magnetic field B = V x A oxt is directed alon g th e 
positive z axis. In 2 + 1 dimensions, there are two inequivalent representations of the Dirac algebra (see e.g.,Ref. [60j): 



7° = cr 3 , 7 1 = iai, 7 2 = ia 2 (D2) 

and 

7° = -er 3 , 7 1 = -*cn,7 2 = -i<?2- (D3) 

which correspond to right- and left-handed coordinate systems. Here Ui are Pauli matrices. The representation of 
gamma matrices (03, icr%, — ia±) used to write the Lagrangian l|2.1|l is related to the representation (03, icri, icr%) by 
means of the unitary transformation U — (I + £<7 3 )/-\/2. Although the final results of the calculations do not depend 
on either the representation of 7-matrices or the gauge, the intermediate expressions depend on this choice. 

Since the representation (|D2|) . (|D3|I is more commonly used in the literature, in this Appendix we solve the Dirac 
equation 1D1I) and obtain the operator (|4.21|) using this representation. 

Let us begin by considering the representation l|D2|) . The energy spectrum in the problem (|D1|I depends on the 
sign of eB; let us for defmiteness assume that eB > 0. Then, the energy spectrum takes the form (to be concrete, we 
assume also that A > 0): 



£ = -M 



(D4) 



E n = ±M n = ±^A 2 + 2n\eB\hv 2 F /c, n=l,2,... 
(the Landau levels). The general solution is 

%/}(x) = 2J KmMnm(l) + b+ m V nm (x)] , (D5) 



where 



■.exp{-iM n t)— ==[ nrj rV-mtt„.m« )' n > 1 ' m<n, 



VOm = — ^= exp(iM t) y j- m ^yim,e J , n = 0, m < 0, (D6) 



Here the functions 61] 



(n + v)\ 



f-i \m 

e -V*Z»/*L»(t), j-™ { $ = A)^ e -W r ™/2 (m < ), (D7) 



where are Laguerre polynomials (L™(£) = for n < —1), I = (Kcj\eB\) x l 2 is the magnetic length, £ = r 2 /2l 2 , 

the quantum number m < n reflects the degeneracy of each level in the angular momentum. The spinors u nrn and 
v nm are normalized as follows: 

d 2 xu ] n , ml (x)u nm (x) = / d 2 xvl i , m ,(x)v nm (x) ^ 5 n ^ n 5 m ' im - (D8) 



Thus the lowest Landau level with n = is special: while at n > 1, there are solutions corresponding to both fermion 
(E n = M n ) and antifermion (E n — —M n ) states, the solution with n = describes only antifermion states. 
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If we used the representation (|D3|I for Dirac's matrices, the general solution would be given by Eq. (|D5|I with 
u nm (x), v nm (x) being substituted by v nm (-x),u nm (-x): 



The solution vo m corresponds to the Landau level n — with positive energy Eq = A, while the solutions v nm , u nm 
with n > 1 correspond to the energy eigenvalues E n = ±M„ (compare with Eq. (D4)). Accordingly, when the spectra 
for two inequivalent representations are united together the degeneracy of the n = level turns out to be a half of 
the degeneracy of the levels with n > 1. The four-component spinor 'J is composed from two two-component spinors 
Eqs. ljD5|l and (|D9I) and using it one can obtain the expression 1)4. 21[) for the pseudospin rotation generator. 
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FIG. 1: (Color online) The optical conductivity a xx (fl) measured in e 2 /h units as a function of the frequency HQ/fj, for two 
different values of B for fi = 50K, T = 15K and T = 10K. We use eB -> (4.5 x I0 4 K 2 )B(T). 
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FIG. 2: (Color online) The Hall conductivity a xy and the diagonal conductivity a xx measured in e 2 /h units as a function of 
the filling v B for T = 2K, T = IK and B = 2T. We use eB -» (4.5 x 10 4 K 2 )B(T) and assume that A = 0. The straight line 
corresponds to a classical dependence a xy = — ecjp|sgn^i/_B [see Eq. IC7I 1. 
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FIG. 3: (Color online) (a) The Hall conductivity a xy and the diagonal conductivity a xx measured in e 2 /h units as a function 
of field B for F = 5K for /i = — 300K and T = 3K. (b) The Hall resistivity p xy = —o- xy /(a xx +a xy ) and the diagonal resistivity 
p xx = o- xx /(a xx + <T xy ) measured in h/e 2 units as a function of field B for T = IK for fi = — 300K and T — 3K. We use 
eB -> (4.5 x 10 4 K 2 )B(T) and assume that A = 0. 
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FIG. 4: (Color online) The Hall conductivity a - ^ and the diagonal conductivity a xx measured in e 2 /h units as a function of 
field _B for F = 6K for fi = 25K and T = 3K. The dash-dotted (black) and dotted (green) lines are calculated using A(B) 
given by Eq. (I4.23H . (b) The Hall resistivity p xy and the diagonal resistivity p xx measured in h/e 2 units as a function of field 
B for F = 6K for /i = 25K and T = 3K. The dash-dotted (black) and dotted (green) lines are calculated using A(B) given by 
Eq. ifOnj. We use eB -> (4.5 x 10 4 K 2 )B(T). 
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FIG. 6: (Color online) The Hall angle Oh as a function of chemical potential /i for two different values of A for B = 3T, 
T = 3K and Y — IK. We use eB -> (4.5 x 10 4 K 2 )B(T). 
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FIG. 7: (Color online) The Nernst signal e y in fiV/K as a function of chemical potential /i for two different values of B for 
T = IK, r = 4K and A = OK. We use eB — > (4.5 x 10 4 K 2 )f?(T). 
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FIG. 8: (Color online) The Nernst signal e y in fiV/K as a function of chemical potential fi for two different values of B for 
T = IK, r = 4K and A = 15K. We use eB -> (4.5 x 10 4 K 2 )B(T). 



